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摘 要 : 本 文 把 Morgado 等 人 提出 的 相互 逼近 (mutual approximation) 技术 从 天 然 卫星 的 天 体位 
置 测量 拓 展 到 小 行星 和 邻近 Gaia 星 表 参 考 星 的 相对 位 置 测 量 ， 和 硕 望 获得 更 高 的 小 行星 的 位 置 测 
量 精 度 。 使 用 中 国 科 学 院 云南 天 文 台 lm 光学 望远镜 在 2020 年 11 月 11~12 日 对 小 行星 702 
(Alauda) 进行 了 试验 观测 。 在 观测 资料 中 ， 目 标 小 行星 共有 5$ 次 相互 逼近 事件 。 参 考 Gaia 
AGE 归 算 结果 显示 ， 在 准确 度 方 面 ， 相 互 逼 近 方 法 与 经 典 天 体 测 量 〈 四 

常数 模型 ) 方法 具有 良好 的 一 致 性 ; 在 测量 精度 方面 ， 相 互 逼 近 方 法 在 赤 经 、 杰 纬 两 个 方向 分 
别 为 6 毫 角 秒 和 2 芝 角 秒 ， 相 比 于 经 典 天 体 测量 〈 四 常数 模型 ) 方法 在 位 置 测量 精度 上 有 6~10 信 
的 提高 。 
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目前 ， 大 部 分 小 行星 的 位 置 测量 都 是 基于 地 面 光学 望远镜 观测 的 图 像 资 料 ， 因 而 提高 地 面 
观测 的 测量 精度 是 一 项 十 分 重要 的 工作 。 人 然而， 影响 小 行星 测量 精度 的 因素 主要 有 两 个 。 一 是 
即使 在 短 曝 光 条 件 下 ， 快 速 运动 的 小 行星 在 CCD 图 像 成 像 上 仍 存 在 拖 影 ， 从 而 引起 测量 精度 的 
损失 。 为 解决 此 问题 ， 文 献 [1] 提 出 了 合成 追踪 技术 (technique of synthetic tracking) 以 获取 准 
确 的 近 地 小 行星 位 置 。 随 后 ，Zhai 等 人 利用 该 技术 在 近 地 小 行星 位 置 测 量 的 精度 上 达到 了 10 毫 
角 秒 说 。 二 是 星 表 提供 恒星 数据 的 精度 问题 。 由 于 星光 在 进入 地 球 大 气 层 后 发 生 抖 动 ， 历 史上 
的 星 表 在 恒星 的 位 置 或 自行 上 精确 度 不 高 ， 这 给 观测 对 象 的 位 置 测量 带 来 了 不 可 忽视 的 误差 。 
随 着 航空 航天 技术 的 进步 与 发 展 ， 人 们 把 探测 器 发 射 到 空间 中 ， 以 消除 地 球 大 气 的 影响 。 欧 洲 
空间 局 成 功 地 发 射 了 Gaia 卫 星 ， 在 完成 阶段 性 测量 任务 后 发 布 了 Gaia 星 表 ， 为 高 精度 的 天 体 测 
量 工作 芮 定 了 重要 基础 。Gaia 星 表 由 于 其 恒星 参考 星 数量 庞大 及 精度 高 的 特性 B494， 在 位 置 归 
算 工 作 中 得 到 国际 同行 的 广泛 认可 。 

2016 年 ，Morgado 等 人 65- 基于 两 天 体 相 对 位 置 测 量 存在 着 精度 回报 (precision premium) [81 
提出 了 相互 逼近 的 天 体 测 量 技术 ， 可 以 精确 地 测量 两 个 天 然 卫 星 的 相对 距离 ， 并 成 功 地 将 该 技 
术 应 用 到 木星 伽利略 卫星 的 天 体 测量 。2019 年 ，Santos-Filho 等 人 也 成 功 地 使 用 该 技术 对 天 王 星 
卫星 进行 了 天 体 测 量 外 。 参 考 JPL (Jet Propulsion Laboratory) 历 表 ， 他 们 的 测量 结果 显示 ， 件 
利 略 卫星 内 部 精度 达到 11.3 毫 角 秒 ， 天 王 星 卫星 达到 45 毫 角 秒 ， 并 且 伯 利 略 卫 星 间 测量 结果 在 
准确 度 和 精度 上 均 与 联合 天 象 技术 相 媲美 。2019 年 ，Lin 等 人 100 基于 M35 星 团 进 一 步 研究 了 精 
度 回 报 现象 ， 发 现 相距 100 角 秒 以 内 的 两 个 天 体 之 间 的 相对 位 置 测量 存在 着 精度 回报 ， 并 且 可 
以 使 用 Sigmoid 函 数 描述 这 两 个 天 体 之 间 的 精度 回报 特性 。 这 项 工作 为 CCD 图 像 上 的 任意 两 个 
天 体 之 间 的 相对 位 置 测量 精度 提供 了 和 定量 描述 。 理 论 上 ， 无 论 两 个 天 体 是 否 发 生 相 对 运动 ， 都 
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存在 着 精度 回报 现象 。 
本 文 把 相互 逼近 技术 从 两 颗 天 然 卫星 天 体位 置 观测 拓展 到 小 行星 和 邻近 恒星 参考 星 〈 位 置 
参考 Gaia 星 表 ) 的 相对 位 置 测量 ， 并 利用 该 技术 以 期 改善 小 行星 的 位 置 测量 精度 。 本 文 的 结构 


如 下 : 第 2 部 分 基于 小 行星 在 较 小 天 空 


星 的 相互 逼近 模型 。 


较 。 最 后 一 部 分 为 本 文 的 总 结 。 


2 相互 逼近 技术 


2.1 理论 模型 
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平面 上 做 色 加 速 直线 运动 的 假设 ， 曾 述 了 小 行星 和 参考 


第 3 部 分 给 出 观测 资料 、 归 算 过 程 和 结果 ， 并 与 经 典 天 体 测 量 方法 进行 比 


如 图 1 所 示 ， 假 设 小 行星 做 匀 加 速 直 线 运动 (加 速度 为 a)〉 经 过 该 天 空 平面 (100 角 秒 以 
经 过 的 路 径 附 近 ， 存 在 一 颗 参 考 星 Ref。 在 观测 时 间 内 ， 小 行星 从 A 点 运动 到 
B 点 ， 并 令 小 行星 在 A 点 的 时 刻 上 = 0《〈 即 起 始点 ) 。 小 行星 到 达 C 点 时 与 参考 星 的 角 距 离 达 


内 ) 。 在 小 行星 


到 最 小 值 《 即 六 


F 击 


(central instant， 用 to 表示 )。 


(例如 图 1 中 的 B 点 ) 。 


显然 ， 角 距离 的 习 


二 中 ， 7 是 小 行星 


=I 


参数 一 一 impact parameter， 用 ao 表示 ) ， 到 达 C 点 的 时 刻 称 为 中 心 时 刻 
对 于 观测 时 间 内 的 任意 时 刻 t， 小 行星 和 参考 星 的 角 距 离 为 dq 


Ref 


区 


1 小 行星 与 参考 星 的 相互 逼近 模型 


Fig.1 Model for the mutual approximation of an asteroid and its reference star 


FEF 方 dq? 有 : 


d’ =d,+(s—s,) =d, +Cvitat /2—s) =a +at+a, t +a +at’ 


\* MERGEFORMAT (1) 


在 时 刻 t= 0 时 的 速度 。 我 们 可 以 使 用 最 小 二 乘法 求解 (1) 式 的 系数 
Cao~Q4 ) , 且 有 关系 : 


\* MERGEFORMAT (2) 


Chinay jy AEH 
Cninax Iva 1 FHH J 


显然 ， 对 (1) 式 求 一 阶 导 并 令 其 等 于 0， 中 心 时 刻 t0 就 是 该 方程 的 根 ， 


a, + 2a,t, +3at +4ati =0 \* MERGEFORMAT (3) 
我 们 可 以 迭代 求解 出 方程 G) 的 根 ， 得 到 中 心 时 刻 t 如 的 值 。 令 上 = 上 一 to， 使 得 中 心 时 刻 t = 0， 
之 后 重新 拟 合 (1) 式 ， 这 时 (1) 式 中 的 参数 so。、to 均 为 0， 于 是 有 关系 : 


40 


qual 


a, = v’, \* MERGEFORMAT (4) 


FH (4) 式 可 知 ， 当 使 用 中 心 时 刻 b = 0 的 观测 数据 再 次 拟 合 〈1) 式 后 ， 可 以 通过 《〈4) 式 直 
接 求解 冲击 参数 do 和 小 行星 在 中 心 时 刻 t0 的 速度 v( 即 小 行星 运动 到 图 1 中 C 点 位 置 的 速度 ) ， 

上 面 的 关系 式 不 仅仅 适用 于 观测 数据 ， 同 样 适 用 于 理论 数据 〈 即 来 自 Gaia DR2 星 表 中 的 
恒星 参考 星 数据 与 历 表 中 的 小 行星 数据 ) 。 因 此 ， 通 过 分 别 拟 合 观测 数据 和 理论 数据 ， 我 们 可 
以 得 到 观测 数据 下 的 中 心 时 刻 女 《单位 为 秒 ) 、 冲 击 参数 d9 (单位 为 像素 ) 和 小 行星 在 中 心 时 
刻 的 速度 v2 〈 单 位 为 像素 / 秒 ) ， 以 及 理论 数据 下 的 中 心 时 刻 妇 《单位 为 秒 ) 、 冲 击 参数 df 
《单位 为 角 秒 ) 和 小 行星 在 中 心 时 刻 的 速度 % (单位 为 角 秒 / 秒 )。 如 果 能 精确 计算 出 CCD 图 
像 的 比例 尺 《〈 见 2.3 节 ) ， 我 们 便 可 以 确定 小 行星 在 中 心 时 刻 的 位 置 偏差 〈O-C， 即 观测 值 减 
去 理论 值 ) : 冲击 参数 的 偏差 ， 即 Ad = pd? 一 d5 (MMAR) ， 其 中 p 是 比例 尺 ， 小 行星 在 
运动 方向 的 位 置 偏差 ， 即 As = vetAt (单位 为 角 秒 )， 其 中 At = (t9 一 区 ) 。 我 们 求解 小 行星 在 
中 心 时 刻 的 速度 虽然 有 两 种 方式 ， 即 来 自 观测 数据 的 pv? 和 来 自理 论 数 据 的 vw， 使 用 任何 一 种 
方式 求解 均 不 影响 观测 资料 的 归 算 结果 ( 见 3.4 市) 。 


2.2 Ad、As 的 投影 


ss 


为 了 将 2.1 节 提 到 的 小 行星 在 中 心 时 刻 的 Ad As 投影 到 赤 经 、 赤 纬 方向 上 ， 我 们 建立 如 
图 2 所 示 的 以 参考 星 Ref 为 原点 的 参考 系 。 图 2 中 ，Ref 表 示 参 考 星 ; &、n 分 别 表示 天 空 平面 
ERZ, WADA; v 为 小 行星 的 运动 方向 ，Ad 和 As 即 为 2.1 节 提 到 的 小 行星 在 中 心 时 刻 


的 位 置 偏差 ; 9 是 小 行星 运动 方向 与 赤 经 方向 的 夹 角 〈( 取 值 区 间 [09,360%])〉。 因 为 Ad 和 As 一 
般 为 不 大 于 0.1 角 秒 的 小 量 ， 我 们 可 以 使 用 历 表 精 确 计 算 9。 将 Ad 和 As 以 向 量 形式 表示 ， 又 人 


7 可 以 视 作 向 量 ( 如 图 2 所 示 ) ， 在 右手 法 则 下 ， 如 果 (Ad x As) AE x ITT AIA, BY BLE 
明 (Ad,As) 和 (Agcos65,A6) 投 影 关 系 为 : 


Aacosé +sin@ cos@ || Ad 
= \* MERGEFORMAT (5) 
Ad —cos@ sin@ || As 


否则 ， 


Aacoso —sin@ cos@ || Ad 
= \* MERGEFORMAT (6) 
Ad +coSO sin@ || As 


图 2 将 Ad 和 As 投影 到 赤 经 、 赤 纬 方向 上 


Fig.2 The projection of the shifts Ad, As to the right ascension and declination 


2.3 比例 尺 


我 们 根据 文献 [11-13] 构 造 CCD 图 像 的 扭曲 模型 ， 去 除 扭曲 效应 后 使 用 四 常数 模型 精确 计 
算 视 场 的 比例 尺 〈 精 度 达到 10-6 角 秒 /像素 ) 。 此 时 由 比例 尺 引 起 pd5 的 误差 对 归 算 结果 的 影响 
其 微 ( 见 3.4 节 ) 。 


3 数据 归 算 与 结果 分 析 


3.1 观测 资料 


表 1 lm 光学 望远镜 和 CCD 参数 说 明 


Table 1 Specifications of the lm optical telescope and CCD chip 


Parameter Value 
Approximate focal length 1330 cm 
F-Ratio 13 
Diameter of primary mirror 100 cm 
CCD field of view 7 xX 7' 
Size of CCD array 2048 x 2048 pixel? 
Adopted scale factor 0.20973 arcsec/pixel 
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2020 € 11 月 11~12 日 ， 使 用 中 国 科 学 院 云南 天 文 台 1 m 光学 望远镜 〈I 滤 光 片 ) 对 小 行 
星 702〈Alauda) 进行 了 观测 ， 观 测 所 用 望远镜 的 参数 说 明 见 表 1。 这 两 个 晚上 一 共 获 得 550 
帧 的 CCD 图 像 数 据 ， 每 天 晚上 拍摄 的 CCD 图 像 数 量 、 曝 光 时 间 和 相互 逼近 事件 的 数量 见 表 2。 
此 外 ， 使 用 I 滤 光 片 进行 观测 可 以 有 效 减 少 较 差 色 差 折射 〈differential color refraction) 效应 对 
天 体 测 量 的 影响 5。 


表 2 小 行星 702 (Alauda) 的 观测 资料 概要 


Table.2 Observations overview for Asteroid 702 Alauda 


: No. mutual 
Date No. of CCD frames Filter Exposure(s) os 
approximation 
2020-11-11 240 I 30 2 
2020-11-12 310 I 30 3 
3.2 数据 归 算 


在 对 参考 星 的 位 置 测量 和 归 算 中 ， 我 们 参考 了 欧洲 空间 局 发 布 的 Gaia DR2 星 表 B4， 小 行 
星 702〈Alauda) 的 理论 位 置 参考 了 美国 国家 航空 航天 局 喷气 推进 实验 室 发 布 的 JPL 历 表 


Chttps://ssd.jpl.nasa.gov) 。 在 上 面 的 位 置 归 算 中 ， 我 们 均 考 虑 到 了 站 心 视 位 置 的 计算 和 大 气 折 
= 射 的 影响 。 
C9 小 行星 与 参考 星相 互 逼 近 的 归 算 步骤 如 下 : 第 一 步 对 所 有 的 观测 资料 ， 使 用 二 维 高 斯 定 心 


算法 拟 合 求解 参考 星 和 小 行星 在 CCD 图 像 中 的 像素 位 置 。 第 二 步 ， 计 算 小 行星 和 参考 星之 间 
象 素 距离 的 平方 ， 并 用 于 拟 合 4 次 多 项 式 函 数 〈 即 (1) 式 ) 。 第 三 步 ， 对 拟 合 得 到 的 (1) 式 
进行 一 阶 求 导 ， 并 令 其 等 于 0《〈 即 方程 G) ) ， 和 迭代 求解 方程 G) 的 根 ， 获 取 中 心 时 刻 。 第 
四 步 ， 将 观测 数据 的 时 间 整 体 减 去 中 心 时 刻 ， 再 次 使 用 观测 数据 拟 合 (1) 式 。 根 据 2.1 节 的 
(4) 式 ， 计算 冲击 参数 和 小 行星 在 中 心 时 刻 的 速度 。 同 理 ， 小 行星 和 参考 星 的 理论 数据 也 进 
行 上 述 类 似 的 归 算 过 程 。 于 是 ， 我 们 便 能 得 到 小 行星 在 中 心 时 刻 的 位 置 偏差 ， 即 2.1 节 的 Ad、 
As。 利 用 JPL 历 表 ， 计 算 小 行星 的 速度 方向 与 杰 经 方 辐 的 来 角 9， 根 据 2.2 节 的 方法 ， 将 Ad 和 
As 投影 到 赤 经 、 赤 纬 方向 ， 最 后 得 到 小 行星 的 位 置 测 量 在 赤 经 、 赤 纬 两 个 方向 的 (0-C) 。 


3.3 结果 与 分 析 


在 2020 年 11 月 11~12 日 拍摄 的 两 个 晚上 的 观测 资料 ， 分 别 存在 2 颗 和 3 颗 参 考 星 与 小 行 
Æ 702 (Alauda) 进行 相互 逼近 ， 其 归 算 结果 见 表 3。 目 标 小 行星 702 (Alauda) 在 JPL 历 表 显 
示 其 星 等 约 为 12mag，5 MES RN G 星 等 范围 为 13~16mag， 于 是 小 行星 和 参考 星 的 星 等 差 
不 大 于 $。 比 较 相 互 逼 近 测 量 技术 与 经 典 天 体 测 量 《〈 四 常数 模型 ) 方法 的 归 算 结果 ( 见 表 3) 
可 以 得 出 ， 两 者 位 置 测量 的 平均 (O-C) 在 赤 经 和 赤 纬 方向 有 良好 的 一 致 性 ， 前 者 的 精度 分 别 
A 6 SAP A 2 毫 角 秒 ， 相 比 于 后 者 有 6-10 倍 的 提高 。 此 外 ， 虽 然 使 用 更 高 阶 模型 的 经 典 天 
体 测 量 方法 可 以 提高 测量 精度 ， 但 是 我 们 在 使 用 相互 逼近 技术 时 ， 并 未 考虑 视 场 的 扭曲 和 高 阶 
大 气 折射 的 模型 ， 因 而 经 典 天 体 测量 (四 常数 模型 ) 方法 与 相互 和 逼近 技术 测量 方法 等 价 。 

图 3 显示 了 表 3 观测 资料 的 归 算 结果 在 赤 经 、 赤 纬 两 个 方向 的 (0O-C) 与 参考 星星 等 和 冲 
击 参数 的 关系 。 从 图 3 可 以 得 出 ， 在 赤 经 、 赤 纬 两 个 方向 的 《0O-C) 与 冲击 参数 没有 明显 的 相 
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Abstract: This paper extends the mutual approximation technique proposed by Morgado et al. from 
positional measurement of natural satellites to a relative positional measurement of an asteroid and its 
nearby Gaia reference star, hoping to obtain a higher positional precision of asteroids. 1m telescope at 
Yunnan Observatory of the Chinese Academy of Sciences is used to obtain the observations of the 
asteroid 702 (Alauda) on November 11 and 12, 2020. There are totally 5 approximation events of the 
asteroid with different reference stars in the observations. With reference to Gaia DR2 (Data Release 2) 
catalogue, our results show that in terms of accuracy, the mutual approximation technique is in good 
agreement with the classical CCD astrometry (i.e., a 4-constant plate model solution). The positional 
precision of the mutual approximation technique is 6 milli-arcseconds and 2 milli-arcseconds in right 
ascension and declination respectively, and the positional precision is improved obviously in 
compared with the classical CCD astrometry by a factor of 6~10. 
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